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The wave- vector (q) and doping (x, y) dependences of the magnetic energy, iron moment, and ef- 
fective exchange interactions in LaFeAsOi-a^Fj, and Bai_2yK2j/Fe2As2 are studied by self-consistent 
LSDA calculations for co-planar spin spirals. For the undoped compounds {x = Q, y = 0), the min- 
imum of the calculated total energy, -E(q), is for q corresponding to stripe antiferromagnetic order. 
Already at low levels of electron doping {x), this minimum becomes flat in LaFeAsOi-xFa; and for 
X > 5%, it shifts to an incommensurate q. In Bai-2i/K2yFe2 As2, stripe order remains stable for hole 
doping up to J/ = 0.3. These results are explained in terms of the band structure. The magnetic 
interactions cannot be accurately described by a simple classical Heisenberg model and the effective 
exchange interactions fitted to -E(q) depend strongly on doping. The doping dependence of the 
iJ(q) curves is compared with that of the noninteracting magnetic susceptibility for which similar 
trends are found. 

PACS numbers: 74.70.-b, 71.20.-b, 75.25.+Z, 75.30.Fv 



I. INTRODUCTION 

The discovery of superconductivity with Tc— 27 K in F- 
doped LaFeAsOi-ajF^; by Hosono and co-workergi a year 
ago initiated an avalanche of experimental and theoret- 
ical investigations of layered iron pnictides and recently 
also chalcogenides. Soon, the superconducting transition 
temperature was raised above 50 K by substituting La 
by smaller rare-earth ionsi^^ The interest in layered iron 
pnictides increased even more when superconductivity 
below Tc—38 K was reported in oxygen-free potassium- 
doped Bai_2yK2yFe2As2r for which good quality single 
crystals could be synthesizedi^ 

Both families of iron pnictides have a quasi-two- 
dimensional (2D) tetragonal crystal structure, in which 
FeAs layers are separated by either LaO or Ba layers. 
The Fe ions form a square lattice sandwiched between 
two As sheets shifted so that each Fe is surrounded by a 
slightly squeezed As tetrahedron. At about 150 K, both 
stoichiometric parent compounds undergo a structural 
transition at which the symmetry of the lattice lowers 
to orthorhombici^ Magnetic order sets in at the same 
temperature as the structural transition in BaFe2As2, 
but at a 20 K lower temperature in LaFeAsO. In both 
cases, the order is striped: ferromagnetic (FM) along the 
shorter axis of the square Fe sublattice and antiferromag- 
netic (AFM) along the longer axis and between the Fe 
layers-i^i^ The Fe moments are 0.4-0.9 in BaFe2As2 
and 0.3-0.4 in LaFeAsO.^-i" 

Electron doping of the FeAs layers in LaFeAsOi_j;Fa; 
suppresses the structural and magnetic transitions in fa- 
vor of superconductivity already at x=0.03r^ Also hole 
doping in Bai_2j,K2j,Fe2As2 suppresses the structural 
and magnetic transition,"* but this requires hole doping 
in excess of y w0.15.*^ Since the superconducting tran- 
sition occurs already for y «0.10, the superconducting 
and striped AFM phases seem to coexist over a fairly 
wide range of hole dopingi^ii^ Although it is not clear 
whether the superconductivity is mediated by AFM fluc- 



tuations or it competes with magnetism, understanding 
merely the magnetic interactions is currently of utmost 
importance. 

A large number of electronic band-structure calcula- 
tions using the local spin-density approximation (LSDA) 
or generalized-gradient approximation now exist for both 
LaFeAsO and BaFcaAss^'^^^'i^^'i^^^^ The re- 
sults obtained for LaFeAsO were reviewed and analyzed 
in Ref. Although the Fe magnetic moment and the 
stabilization energies of different magnetic solutions de- 
pend strongly on the computational method and the 
exchange-correlation functional, as well as on whether 
the experimental or calculated structure is used, all cal- 
culations predict that stripe AFM order is the mag- 
netic ground state in both parent compounds. However, 
the calculated sublattice magnetizations are significantly 
larger than the ones deduced from the neutron diffrac- 
tion, /iSR, and Mossbauer experimentSi^i^iSiii^ 

For both parent compounds, many authors (see, e.g., 
Refs. 11,„,16;,20,) have noticed a strong Fermi surface (FS) 
nesting for the q vector which corresponds to stripe AFM 
order between the Fe d^z/yz'^^^^ hole sheet and one of the 
two electron sheets. This nesting causes peaks in both 
the imaginary and real parts of the non-interacting spin 
susceptibility, Xoiq), at the stripe q. Although electron 
doping of LaFeAsO suppresses the peak and shifts it to 
an incommensurate wave vectorj^ it is widely believed 
that stripe AFM order remains the LSDA ground state 
of LaFeAsOi_2;Fa;. That merely 3% electron doping suf- 
fices to destroy the static stripe order, has been related to 
filling of the three-dimensional 3D Fe d^z^-i-Mke band.— 
However, calculations which use the experimental struc- 
ture -such as those presented below- place the top of the 
f^Sz^-i-like band several hundred meV below the Fermi 
level. 

In this paper, we shall present LSDA calculations of 
moments and energies of spin spirals in LaFeAsOi_a;Fa; 
and Bai_2i/K2iyFe2As2 as functions of wave vector, q, and 
doping, X or y. We find that upon increasing doping. 
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TABLE I: Lattice and LMTO spheres 



Sphere Wycoff position 


Radius (A) 


LaFeAsO {P4/nmm), a^b = 


4.04 A, c = 8.74 A 


Fe 2b (0.25, 0.75, 0.50) 


1.41 


As 2c (0.25, 0.25, 0.651) 


1.48 


La 2c (0.25, 0.25, 0.142) 


1.71 


2a (0.25, 0.75, 0) 


1.13 


El 2c (0.25, 0.25, 0.407) 


1.08 


E2 4/ (0.25, 0.75, 0.266) 


1.04 


E3 2c (0.25, 0.25, -.115) 


0.97 


BaFe2As2 (M/mmm), a = b = 


3.96 A, c = 13.02 A 


Fe 4d (0, 0.5, 0.25) 


1.38 


As 4e (0, 0, 0.355) 


1.46 


Ba 2a (0, 0, 0) 


2.01 


El 4e (0, 0, 0.195) 


0.99 


E2 2b (0, 0, 0.50) 


0.78 


E3 8g (0.5, 0, -.105) 


0.84 



stripe AFM order becomes unstable in favor of an in- 
commensurate spin-density wave (SDW). 

Before getting to the spin spirals, we shall explain 
our computational method and compare with previous 
results -and also present results involving spin-orbit 
coupling- for the paramagnetic band structures and the 
commensurate stripe and checkerboard SDWs for the 
parent compounds. 

II. COMPUTATIONAL METHOD 

Our LSDA scalar-relativistic calculations for co-planar 
spin spirals in LaFeAsOi_a;Fa; with x=0, 0.1, 0.2, and 
0.3 and Bai_2yK2yFe2As2 with y—Q, 0.1, 0.2, and 
0.3 were carried out using the linear muffin-tin or- 
bital (LMTO) method in the atomic-sphere approxi- 
mation and including the combined correction termi^ 
For the exchange-correlation potential, we used the 
Perdew-Wang parametrization.^'* Charge- and spin-self- 
consistent calculations were carried out for all doping 
levels using the virtual-crystal approximation with frac- 
tional atomic number of O or Ba. All our calcula- 
tions were for the experimental room-temperature crys- 
tal structures of the undoped compounds ^^i^ This is im- 
portant. Whereas the primitive cell of LaFeAsO holds 
two formula units and is tetragonal (PA/nmm), that of 
BaFe2As2 holds one unit and is body-centered tetrago- 
nal (bet) {lA/mmm). The lattice constants as well as the 
positions and radii of the space-filling atomic and empty 
(E) LMTO spheres are given in Table [J 

For the Fe d orbitals, such as dxy^ we use the x and 
y axes (not to be confused with the levels of electron 
and hole doping) which span the quadratic Fe sublattice 
and are therefore turned 45° with respect to the tetrag- 
onal a and h axes. The lobes of the Fe dxy orbital thus 
point toward the projections of the As sublattice onto 
the Fe plane, while dj.2_y2 points toward the nearest Fe 
neighbors. Or in other words: the Fe-Fe ddir interaction 




FIG. 1: (Color online) FS cross sections (fcz=0) for LaFeAsO 
in the tetragonal Brillouin zone (BZ). The a direction is hor- 
izontal, the b direction is vertical, and the coordinates of the 
(black) F, X, and M points are as given in Fig. [1 The FS 
on the left- and right-hand sides were calculated, respectively, 
with and without a 150 meV downwards adjustment of evxy 
in the LMTO method. The F-centered d^y-like hole cylin- 
der is shown in thin (red) line and the concentric dj^^^zx-like 
hole cylinders are shown in thick (blue) lines. The two M- 
centered, electron cylinders shown in thin (green) lines have 
super-ellipsoidal cross sections with main axes directed to- 
ward r. The fact that a primitive translation of the square 
Fe sublattice followed by mirroring in the Fe plane generates 
an Abelian subgroup of the space group allows one to fold the 
band structure out to a large one-formula-unit BZ. This moves 
the d^z/yz-like hole cylinders to the nearest-neighbor P points, 
which are the (blue) M corners of the large BZ and have the 
coordinates (1, l)7r\/2/a in the (a;, 2/)-system. This also sepa- 
rates the electron cylinders onto different M points, now called 
X (1,0) and Y (0,1) (blue), such that the super-ellipses have 
their short axis pointing toward the F-centered d^^y-like hole 
cylinders and the long axis pointing toward the M-centered 
c^rrz/j/z-like hole cylinders. Introducing a commensurate SDW 
with q = (1,0) in the {x,y) system, i.e., FM-ordered stripes 
in the j/-direction and AFM order in the a::-direction, will shift 
the hole cylinders by q. This is shown by dashed lines around 
the bottom-left and top-right corners, M=X. The selection 
rules following from the glide mirror allow coupling only be- 
tween the dxy holes and the electron sheet with the short 
axis along q, and between the d^^/y^ holes and the electron 
sheet with the long axis along q. Introducing a SDW with 
q = (1,1), i.e., checkerboard AFM, folds the large BZ back 
into the small tetragonal one. Here the nesting is between the 
da^y-like and one of the d^-z/yz-like hole sheets and between 
two different electron sheets. Not included in our spin-spiral 
calculations is the spin-orbit coupling (^ ~ 60 meV), which 
invalidates the glide mirror and violates the selection rules. 
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FIG. 2: (Color online) Band structures of undoped LaFeAsO (left panel) and BaFe2As2 (right panel) obtained from non-spin- 
polarized calculations for the experimental structure. The size of (red) circles and (blue) squares is proportional to, respectively, 
the Fe d^y and dy^^zx partial- wave character. The coordinates shown at the bottom are in the tetragonal (a, b, c)-system where 
the cell contains two Fe atoms, i.e., in units of (1/a, l/6=l/a, 1/c). The Fermi level is the zero of energy. 



involves the dxy orbitals and the Fe-Fe dda interaction in- 
volves the dx2_y2 orbitals. This convention is the same as 
the one used for the Cu orbitals in the high-temperature 
superconducting cuprates. The distance, a/\/2, between 
Fe nearest neighbors is 2.85 A in LaFeAsO and 2.80 A in 
BaFe2As2. 

Our band structures and Fermi surfaces for undoped 
paramagnetic LaFeAsO and BaFe2As2 obtained from 
spin-restricted LDA LMTO calculations agree well with 
those obtained with the same method and presented in 
Ref. [m. 

However, for LaFeAsO the LMTO Fermi surface dif- 
fers from that obtained with the full-potential linear aug- 
mented plane-wave (LAPW) method;^ Whereas LMTO 
finds the two innermost F-centered hole cylinders to be 
dxz,yz hke and the outermost dxy hkc, LAPW finds the 
opposite order. Since LAPW is computationally more 
cumbersome and accurate than LMTO, we used LMTO 
but applied an external crystal field which shifted the en- 
ergy (e^xy) of the Fe dxy partial wave downwards by 150 
meV. This adjustment brought the LMTO and LAPW 
band structures into almost complete agreement. The 
adjusted and unadjusted Fermi surfaces are shown re- 
spectively on the left and right-hand sides of Fig. [T] In 
the following, all results presented for LaFeAsO i_a;Fj: ~ 
such as the paramagnetic bands on the left-hand side of 
Fig.[2]and the spin-spiral moments and energies in Fig.[4]- 
are those obtained with LMTO and the dxy-energy down- 
shifted, unless otherwise stated. 

For BaFe2As2, the LMTO LDA band structure shown 
on the right-hand side of Fig. [2] agrees very well with 
the LAPW one. For that reason, our LMTO calcula- 
tions for Bai_2j,K2yFe2As2 were all performed without 
any adjustment. 



Finally, we calculated effects of spin-orbit cou- 
pling (Fermi surface splittings and magnetocrystalline 
energies) using the spin-polarized relativistic LMTO 
method-- and the experimental structure. 



III. PARAMAGNETIC ENERGY BANDS 

The paramagnetic scalar-relativistic bands for un- 
doped LaFeAsO are shown near the Fermi level at the 
left-hand side of Fig. ^ We see that at F, the top of the 
dxz/yz band (blue squares) is ~180 meV while that of 
the dxy band (red spheres) is merely ~30 meV above ep- 
Without the 150 meV shift, these levels would have been 
nearly degenerate, and since the mass of the dxy band is 
higher that those of the two dxz/yz bands, the cylindrical 
dxy sheet would have been the widest, as shown in the 
right-hand side of Fig. [T] Although spin-orbit coupling 
cannot be included in spin-spiral calculations, we men- 
tion that it splits the degenerate top of the dxz/yz-^^^^ 
band by 50 meV, a value consistent with ^PoSd ~ 60 meV 
and a 20% reduction due to by-mixing of As p character. 

With an even number of electrons (Fe d^), the sum of 
the volumes of the three hole sheets equals that of the 
electron sheets. Of those, there are two equivalent cylin- 
ders centered at M and turned 90° with respect to each 
other (thin green). Their cross sections are superellip- 
soidal with main axes pointing toward F. The superellip- 
soidal cross section arises because these sheets result from 
a dxz or dyz band hybridizing with a lower-lying dxy-^ike 
band.~ The main character is dxz or dyz near the short 
axis and dxy near the long axis. These electron cylinders 
have more As p character and more fcz-dispersion than 
the hole cylinders, and are therefore more warped. 
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The fact that a primitive translation of the square Fe 
sublattice followed by mirroring in the Fe plane generates 
an Abelian subgroup of the Pi/nmm space group allows 
one to fold the band structure out to the large Brillouin 
zone (BZ) well known from the cuprates.^^'^^ This folding 
out has the advantage of separating the dxy cylinders 
from the d^z/yz ones by placing them at respectively F 
(0,0) and M (1,1) in the [x^y) system. It also separates 
the two electron pockets from each other by placing them 
at, respectively, X (1,0) and Y (0,1) with the long axis 
pointing toward M. 

Note that in order to distinguish the tetrago- 
nal (a, b, c)-directed coordinates in reciprocal space 
from the quadratic (a;, y, z || c)-directed ones, we use 
the units (1/a, 1/a, 1/c) in the former case and 
(-\/27r/a, -\/27r/a, tt/c) in the latter case. Hence, tetrag- 
onal (quadratic) reciprocal-space coordinates are recog- 
nized by the presence (absence) of the factor tt. 

Spin-orbit coupling does not commute with the above- 
mentioned glide mirror and will therefore split the cross- 
ing between the dxy and one of the d^z/yz hole bands 
by about 50 meV, as well as the crossings between two 
electron bands as has been observed in LaFePOi^ 

At the right-hand side of Fig. [51 we show the paramag- 
netic bands near the Fermi level for undoped BaFe2As2. 
Since the As-Fe2-As layers are separated by a thin Ba 
layer, rather than by a thick La-02-La layer, and the AS 
atoms along the z axis are stacked on top of each other 
the band structure of BaFe2As2 disperses more in the z 
direction than that of LaFeAsO. But apart from that, the 
band structures are very similar. For ease of comparison 
with the LaFeAsO band structure at the left-hand side 
of the figure, we have chosen the same route in (fc^,, ky) 
space, but have taken — / (2c) -except in the very 
last panel" because this choice makes the amplitude of 
the Bloch waves vanish in the Ba-plane and thus min- 
imizes the effects of kz dispersion. Along (0,0, /cz), the 
(ia;j,-like band is 140 meV above e f and does not disperse 
with kz^ whereas the doubly degenerate d^z/yz-^^^^ band 
disperses from 40 to 200 meV above ep with kz increas- 
ing from to tt/c, and thus goes from below to above 
the top of the d^y band. The hole sheets thus remain 
cylinders although the d^z/yz'^^^^ sheets are significantly 
warped. Near (0,0, tt/c), we see a band with strong kz- 
dispersion dip below ep- This band is As Pz like and 
cannot hybridize with the d^y-hke band, but only with 
the dxz/yz-^^^^ band, but not along (0,0, fcz). This hy- 
bridization gives rise to an intricate shape of the d^z/vz 
bands near e^? for fcz 7^ tt/ (2c) and is discussed in Ref.j2Q. 
A further difference with the LaFeAsO bands is that the 
electron cylinders around (7r,7r,fcz) are not degenerate 
along the a and b directions and that the kz dispersion 
of the dxy-\ike component is as large as 150 meV due to 
by-mixing of As Pz character. 

The main effect of electron doping on the band struc- 
ture is to move the Fermi level up or -equivalently- to 
move the bands down with respect to the Fermi level. 
This is clearly seen in Fig. [3] where we show dependence 
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FIG. 3: (Color online) Band positions as functions of doping 
in the range —30% < x = —y < 30% for Bai_2yK2i,Fe2As2 
(left) and LaFeAsOi-^cF^ (right). 



on doping - ranging from —30% in Bai_2iyK2i/Fe2As2 to 
-f 30% in LaFeAsOi_a;Fa; - of the top of the da;j,-like (red 
points) and dj-^/jj^-like (blue squares) hole bands at (0, 
0, 7r/2, or 0) and the bottoms of the d^y- and dxz- or 
dyz-like electron bands at (tt ,7r ,7r/2, or 0). We note, 
first of all, that the dominating rigid shift of these levels 
is roughly continuous when passing from hole doping in 
Bai_2yK2iyFe2As2 to electron doping in LaFeAsOi_a;Fa; 
(and even more so had we not corrected e^, xy for the lat- 
ter). Second, we note that 11% electron doping fills the 
dxy band, and 33% fills the dxz/yz band, and that no 
other Fe c?-like band gets filled or emptied in the ±30% 
doping range. Within this range, the bands move by 
about 300 meV, corresponding to an average density of 
states of order 2 electrons/Fe/eV. The deviation from 
rigid-band behavior, i.e., parallel movement of all levels, 
is roughly 50% of this. 



IV. MAGNETIC MOMENTS AND ENERGIES 
OF COLLINEAR, COMMENSURATE SDWS IN 
THE PARENT COMPOUNDS 

For LaFeAsO, our LMTO calculations with imposed 
FM order converge to a nonmagnetic solution. For stripe 
and checkerboard AFM orders we calculate magnetic sta- 
bilization energies of, respectively, 78 and 39 meV/Fe, 
values which compare well with those, 84 and 26 meV/Fe, 
obtained from LAPW calculations for the experimental 
structure. Our sublattice magnetizations (Fe moments) 
of 1.3 for stripe and 1.2 /ie/Fe for checkerboard AFM or- 
der are somewhat smaller than those, 1.8 and 1.5 ^e/Fe, 
found with LAPW.^^ This could be due to integrating 
the spin-density over different regions, i.e., LMTO and 
LAPW spheres have different sizes. The LSDA band 
structure and Fermi surface for stripe order are in a 
good agreement with those obtained with LAPW.— Fi- 
nally, we mention that without the downwards e^xy shift. 
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LMTO yields slightly larger energies, 95 and 53 meV/Fe, 
and moments, 1.4 and 1.4 ^e/Fe. 

For BaFe2As2, we also find no FM solution. For stripe 
order, we find a somewhat lower stabilization energy, 62 
meV/Fe, than for LaFeAsO but the same moment, 1.3 
Aie/Fe. Using LAPW (Ref Hi) instead of LMTO, we 
find the same stabilization energy, 85 meV/Fe, as for 
LaFeAsO, and a slightly smaller moment, 1.7 /ie/Fe. 
The calculated LMTO and LAPW band structures for 
the striped phase agree reasonably well. LMTO calcu- 
lations performed for the experimental low-temperature 
orthorhombic structure^, with the FM stripes oriented ei- 
ther along the longer or the shorter axis result in a lower 
total energy for the latter. This is in accord with the 
experimental data in Ref. and a previous calculation;— 

For stripes in both LaFeAsO and BaFe2As2, the max- 
imum exchange splitting, i.e., that of degenerate bands, 
Ejk and eyk+q, with the same dominant d character 
is A « 1.2 eV, and this is consistent with an Fe moment 
of M w 1.3 /xb/Fc and Stoner theory: A = MI, with 
the usual value of the Fe effective Stoner parameter (~ 
Hunds rule Jh)- I « 0.9 eVi^i^ Since A is an order 
of magnitude larger than the energies of the (paramag- 
netic, un-coupled) electron and hole pockets with respect 
to the Fermi energy, but considerably smaller that the 
sub-band widths -which are several eV- the strength of 
the LSDA exchange coupling is intermediate when the 
magnetic order is stripe or checkerboard. This means 
that neither Fermi surface nesting, nor the difference be- 
tween the band structures of pure LaFeAsO and pure 
BaFe2As2, nor the 150 meV adjustment for LaFeAsO, 
strongly affects the magnetic moments and energies. 

Although the spin-orbit coupling in the Fe 3d shell is 
weak (^Fo3(i ~ 60 meV), it does lead to a distinct mag- 
netocrystalline anisotropy with the Fe moment lying in 
the plane and perpendicular to the FM stripes. This 
has been shown by neutron diffraction for BaFe2As2i- 
Our relativistic LMTO calculations agree with this: For 
LaFeAsO we find that it costs 0.27 meV/Fe to turn the 
moment in the plane from the easy to the FM-stripe di- 
rection, and an additional 0.13 meV/Fe to turn the mo- 
ment perpendicular to the plane. For BaFe2As2, the cor- 
responding energies are 0.16 and 0.04 meV/Fe, i.e., the 
anisotropies are smaller. 



V. MAGNETIC MOMENTS AND ENERGIES 
OF SPIN SPIRALS IN THE DOPED 
COMPOUNDS 

A. Spin spirals 

A spin spiral is characterized by the following proper- 
ties: Upon a lattice translation t, the magnitude of the 
magnetization density and its projection onto a global z 
direction are unchanged, but the projection onto the per- 
pendicular (x, y) plane rotates by an angle ip (t) propor- 
tional to the projection of the translation onto the wave 



vector q of the spin spiral. The magnetization-density of 
a spin spiral thus satisfies the equation 

M (r + t, q) = M (r, q) 

• {zz cos 9 + [xx cos f{t) -I- yy sin 1^ (t)] sin 9} 
with (f (t) = Trq • t. 

(The factor tt merely follows from the convention chosen 
in Sec. Illll for t referring to the Fe sublattice.) Examples 
of such spin spirals are shown by the solid arrows in Fig. [4] 
at the bottom of the right-hand side. Note that for Fe 
atoms lying along rows perpendicular to q, the alignment 
is FM. 

In order to solve the one-electron problem in the pres- 
ence of such a spin spiral, one may span the one-electron 
Hilbert space by a complete set of localized orbitals, 
(/)j(r — t), times pure spin-functions, Xt (f) = at (c) 
or Pt (f ) , whose quantization direction is chosen along 
the local direction of the magnetization. In this rep- 
resentation, the one-electron Hamiltonian without spin- 
orbit coupling is translationally invariant, albeit with q- 
dependent hopping integrals, so that there is no coupling 
between Bloch sums, 0j (r — t) xt (f) exp (Trik • t), 
with different Bloch vectors. As a consequence, the one- 
electron problem can be solved for any q, without increas- 
ing the size of the primitive cell, provided that spin-orbit 
coupling is neg lected,^^ This, together with the LSDA, 
enables simple calculation of spin-spiral moments and to- 
tal energies. 

In our LMTO calculations, the localized orbitals were 
taken to be the partial waves truncated outside their 
sphere. This means that we forced the direction of mag- 
netization to be constant inside each sphere. The mo- 
ment that we quote is the one integrated over an Fe 
sphere. We considered spin spirals for which the mag- 
netization is in the Fe {x,y) plane, i.e., 9 = tt/2 and first 
took q to lie in the plane and t to span the square Fe 
sublattice. To achieve this in calculations employing the 
tetragonal translational cell, the phases ipi, which deter- 
mine the magnetization directions in two Fe spheres at 
positions , were fixed to (pi — q - ti. 

The spin spiral with q = f (0, 0) produces FM order. In 
the spin spiral with q = X(l, 0), the moments rotates by 
TT upon translation by t =(1,0), and by upon transla- 
tion by t = (0, 1), i.e., the order is stripe with AFM align- 
ment of nearest-neighbor moments along the x direction 
and FM alignment along the y direction. In the spin 
spiral with q = M (1, 1), the moments rotate by tt upon 
translation by t = (1, 0), as well as by t = (0, 1), i.e., the 
order is checkerboard with all four nearest-neighbor mo- 
ments antiparallel and all four second-nearest moments 
parallel. These spin spirals with q at high-symmetry 
points are all coUinear and commensurate. 

NoncoUinear and incommensurate spin spirals with q 
near -but not at- X and M are illustrated at the bottom 
of the right-hand side of Fig. |4l When q is on the FX 
line, the order along the y direction remains FM, while 
going from one Fe to the next in the positive x direc- 
tion, the moment rotates by between and tt. When 
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FIG. 4: (Color online) Magnetic moments (upper panels) and energies (lower panels) per Fe of in-plane spin spirals as functions 
of q = {qx, Qy, 0) along the boundaries f (0, 0) — X (1, 0) — M (1, 1) — F (0, 0) of the irreducible part of the large BZ for different 
electron (x) and hole (y) dopings. The energies of the ji-j2 Heisenberg model for x{y)=0 and 0.2 are given by, respectively, 
dashed and dash-dotted lines. Representative real-space spin structures are shown at the bottom right for the q vectors denoted 
by dots. 



q is on the XM line, the order along the x direction is 
AFM, while upon going from one Fe to the next in the 
positive y direction, the moment rotates by between 
and TT. When finally q is on the MF line, the moment 
rotates by the same angle, lying between and tt, regard- 
less of whether t — (1,0) or (0, 1), and the order along 
the second-nearest-neighbor direction perpendicular to q, 
i.e., for t = (1,-1), is FM. 

We now discuss the calculated results shown in this 
figure for the magnetic moments (top) and energies (bot- 
tom) of spin spirals in electron- (left) and hole-doped 
(right) compounds as functions of q along the triangular 
boundary F-X-M-F of the irreducible part of the large 
BZ (see Fig. [J). 



B. Pure and electron-doped LaFeAsOi_i,F^ 

For undoped LaFeAsO, i.e. for x = (red dots), the 
lowest energy is reached at the X point, i.e., for stripe 
AFM order. This agrees with the results of previous 
calculationsi^ii^i2£ and experimental dataJ^ When mov- 
ing from X(l, 0) toward f (0, 0), the AFM order between 
nearest-neighbor Fe rows along y is destroyed, and this 



leads to a rapid decrease in moment and increase in en- 
ergy; for |q| = g < 0.6 the self-consistent solution is non- 
magnetic. When moving from X(1,0) toward M(l, 1), 
the FM order in the Fe rows along y is destroyed and be- 
comes AFM at M. Whereas the moment at first remains 
nearly constant at 1.4 /ie, but finally decreases to 1.2 
fiB at M, the energy increases nearly parabolically from 
-78 to -39 meV. Moving from M(l, 1) toward f (0, 0), 
the nearest-neighbor AFM order develops toward FM or- 
der by preserving the FM order between second-nearest 
neighbors along t = (1,-1) but destroying the one along 
t = (1, 1). Hereby the moment first increases slightly, but 
then decreases rapidly and vanishes when q < 0.7. The 
energy falls to a local minimum at q « 1, and then in- 
creases rapidly for q decreasing to 0.7 where the moment 
disappears. 

As shown by the dashed curve in Fig.|4l for q > 0.7 our 
calculated E (q) for undoped LaFeAsO is approximated 
reasonably well by a classical Heisenberg model on the 
square lattice with AFM exchange coupling constants ji 
and j2 between, respectively, nearest and next-nearest 
neighbors. We obtain the ratio j2/ji = 0.71 (> 1/2) by 
fitting to the q position of the local minimum along MF 
and then obtain the values ji = 81 and j2 = 57 meV 
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by fitting the energy difference i?min (MF) — E (X) and 
the value S = (X) = 0.67. The q-independent con- 
stant is finally chosen such that the Heisenberg model 
fits the calculated energy at X (and the minimum along 
Mr). Our values of ji and j2 are comparable to those 
obtained from Fig. 3 in Ref. [3 when interpolated to 
M (X) = 1.3 /iB, although our j^/ji is somewhat higher. 
The anisotropic exchange coupling constants calculated 
in Ref. [l^ are meaningful only for small deviations from 
stripe AFM order, that is for q near X, and cannot be 
directly compared with our effective ji and j2- Neverthe- 
less, our jiS"^ = 420 K and j^-S^ = 297 K are of the same 
order of magnitude as —J^ — 550 K and — = 260 
K of Ref. [ij. This indicates that the three different 
approaches result in the comparable exchange interac- 
tions. Although the overall shape of the dashed E (q) 
curve given by the Heisenberg model is similar to the 
red-dotted one obtained from our LSDA calculation, dis- 
crepancies can be clearly seen even in the part of the BZ 
where the Fe moment remains nearly constant: The cal- 
culated energy is higher in the vicinity of the M point and 
it increases far more rapidly when going from X toward 
f. 

Our calculations reveal that a new local minimum of 
E{q) develops along XM near (1,0.3) when the elec- 
tron doping exceeds about 5%. For 0.10 < x < 0.25, 
this minimum is deeper than the one a X, i.e. it is 
the global minimum, and this means that stripe order 
is unstable against formation of an incommensurate non- 
collinear SDW. This is in line with experimental phase di- 
agrams for REFeAsOi-ajFj; compounds which show that 
magnetic order is rapidly suppressed by F doping ;ii For 
X — 0.3 we find that the energy minimum has returned 
to X, although a well-defined shoulder can still be seen 
at (1,0.3). This instability toward an incommensurate 
SDW cannot be reproduced by fitting to the ji-j2 Heisen- 
berg model, as is clearly seen by comparison of the dot- 
dashed and blue T curves. 

Electron doping is seen to increase the moment for 
stripe AFM order from 1.4 to 1.7 for x = 0.3. But 
this increase is localized to q being close to X. 

The destabilization of stripe AFM order with electron 
doping seems to be caused by occupation of a narrow 
peak of the density of states (DOS), which in undoped 
LaFeAsO is ^150 meV above the Fermi level (ep)- The 
band responsible for this peak is the paramagnetic Fe dyz 
band which hardly disperses in the direction and stays 
within ±200 meV of the Fermi level over a region of k 
space near the entire MY line. This is the band seen in 
Fig- m to connect the dy^-like hole pocket at M with the 
superellipsoidal electron pocket at Y. Now, FM stripes in 
the y direction with AFM order along x will couple states 
at k with those atk-|-q = k-|-(l,0), i.e. will fold the 
large BZ perpendicular to the y direction, thus placing M 
on top of Y, and exchange-split states with the same Fe 
d character by ±0.5 eV. As a result, the flat dyz band 
will have its upper minority-spin dyz band ~150 meV 
above in the undoped compound. Since the DOS 



near ep is very low in the stripe-ordered undoped com- 
pound, not much electron doping is needed to occupy 
part of the flat dyz minority-spin band. For q moving 
away from X in the perpendicular direction, i.e., toward 
M, the DOS peak soon splits in two. The concomitant 
decrease in band energy for dopings such that sp is in 
the valley between the subpeaks compensates for the de- 
crease in negative magnetic exchange energy (— ^M^/) 
caused by the decrease in magnetization seen in the up- 
per left-hand side of Fig. [D The magnetization decreases 
because the nesting is less good (less phase space avail- 
able for gapping) when q moves so far away from X in 
the perpendicular direction that q does not so well trans- 
late the k tube around MY in which the dyz band is flat 
onto itself. The energy minimum finally shifts back to 
q = X (1, 0) once the electron doping is so high (> 25%) 
that the narrow DOS peak is completely filled, i.e., when 
Ei? is above the flat part of the minority-spin dyz band. 

Recently, a commensurate-to-incommensurate transi- 
tion with increasing x in LaFeAsOi^^^F^; has been re- 
ported in Ref. fssi on the base of full-potential LAPW 
spin-spiral calculations. These calculations were however 
performed using calculated As positions and can there- 
fore not be compared directly our calculations based on 
the experimental structure. 

Returning now to the undoped compound and q mov- 
ing away from X toward M, sp is below both subpeaks 
and the Heisenberg-type dependence of E{q) persists un- 
til in the vicinity of M a huge DOS peak appears just 
below ep. 

For stripe-ordered pure LaFeAsO, another narrow 
DOS peak exists 270 meV above ep and arises from the 
upper minority-spin c?2;j,-like band being split from the 
lower minority-spin band due to folding of the M-centered 
hole pocket onto the Y-centered electron pocket;2£ 

We estimated the strength of the interlayer exchange 
coupling by performing calculations for spin spirals with 
q not lying in the Fe plane. Specifically we took q = 
{1,0, Qz) corresponding to stripe order in each layer and 
rotation of the magnetic moments between layers. These 
calculations resulted in a very small change in energy 
when changing the alignment of the Fe moments along 
the c direction from ferromagnetic to antiferromagnetic. 
A weak dependence of the energy on the magnetic order 
along the c axis was also reported in Ref. [l^. This nearly 
two-dimensional character of the magnetic interactions is 
not affected by F doping. 

The q dependencies of the spin-spiral moment and en- 
ergy calculated for LaFeAsOi_a;Fj; without downshifting 
the dxy energies are shown in Fig. O Also in this case is 
the energy minimum of undoped LaFeAsO at X, but ^ 20 
meV deeper. The local minimum along ML is at the same 
q position, but merely 10 meV deeper, so that the fitted 
values, ji — 93 and j2 = 66 meV, are a bit larger but 
have the same ratio. For electron-doped LaFeAsOi_a;Fj;, 
the q dependence of the energy is somewhat stronger and 
exhibits a shoulder at q w (1, 0.3) which does not develop 
into a well-defined minimum. Nevertheless, destabiliza- 
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FIG. 5: (Color online) Same as Fig.|4l but for LaFeAsOi-iFa; 
without downshifting the Fe dxi/ energies. 



TABLE II: The doping dependence of the Fe magnetic 
moments A4"(X) calculated for stripe AF and the ex- 
change coupling constants ji and in LaFeAsOi-^Fa; and 



Bai_2aK2i;Fe2As2. 


For LaFeAsOi 


^xFx, the values calculated 


both with and without 


the shift 


of the dxy 


states 


are pre- 


sented. The latter 


are g 


iven in parentheses. 






Compound 


x/y M(X) (MB 


) h/ji ji 


(meV) j2 (meV 


LaFeAsOi-j;Fx 


0.3 


1.68 


0.71 


61 


43 






(1.74) 


(0.71) 


(77) 


(54) 




0.2 


1.58 


0.71 


47 


33 






(1.65) 


(0.62) 


(103) 


(63) 




0.1 


1.47 


0.66 


71 


47 






(1.56) 


(0.59) 


(160) 


(94) 







1.34 


0.71 


81 


57 






(1.42) 


(0.71) 


(93) 


(66) 


Bai-2j;K2j,Fe2As2 





1.34 


0.77 


95 


73 




0.1 


1.19 


0.85 


77 


65 




0.2 


1.11 


1.10 


39 


43 




0.3 


0.99 


1.10 


19 


21 



tion of the commensurate collinear stripe order does oc- 
cur, although the details are seen to depend sensitively 
on the underlying band structure. 



C. Pure and hole-doped Bai 2yK2j;Fe2As2 

The spin-spiral moments and energies as functions of 
q = {qx,qy,0) for Bai^2yK2yFe2As2 with ?;=0, 0.1, 0.2, 
and 0.3 holes per FeAs are shown in the right-hand 
panel of Fig. [H For undoped BaFe2As2, these curves 
(red dots) are qualitatively similar to those calculated 
for LaFeAsO, and the trend for increasing hole-doping 
of Bai_2j;K2j,Fe2As2 continues the trend for decreasing 
electron-doping of LaFeAsOi_a;Fa;, i.e., y ^ — x. The cal- 
culated energy gain E(T) — E{X.) due to the formation of 
stripe AFM order in BaFe2As2 is somewhat smaller than 
in LaFeAsO, but the energy difference E^i^ {Mr)-E (X) 
is a bit larger and this causes the fitted values ji = 95 
and j2 — 73 meV to be a bit larger. The ratio j2/ji=0.77 
is slightly larger than for LaFeAsO. Hole doping strongly 
reduces the stripe- formation energy, E{r) — E{X). Nev- 
ertheless, the energy minimum remains at X until the 
hole doping exceeds 25%, at which point the minimum 
splits in two with the lowest lying along XF. Our calcu- 
lations thus show that stripe order is more resistant to 
hole doping in Bai_2i/K2yFe2As2 than to electron doping 
in LaFcAsOi_a;Fa;. This conclusion is supported by ex- 
perimental observations of the traces of the spin-density- 
wave phase for K-doping as high as 0.4, i.e., well into the 
superconducting region^ii^ 

The local minimum along MF moves toward F as hole 
doping increases. As a consequence, the estimated j2/ji 
increases and reaches the value of 1.1 when y = 0.2. Since 
the energy difference i^inin (MF ) — E (X) decreases with 
hole doping, the values of the effective coupling constants 
decrease to ji — 39 and 43 meV when y = 0.2. 

The Fe moment calculated for stripe AFM order de- 
creases from 1.34 /iB in the undoped compound to 0.99 
fiB for y = 0.3. The maximum of M(q) is however not 
at X, but at q « (1,0.4) along XM. This maximum 
becomes more pronounced with hole doping. In con- 
trast to the situation in LaFeAsOi-a^Fj;, where for q in 
a large region around F the non-magnetic solution is sta- 
ble, in Bai_2yK2j,Fe2As2 magnetic solutions exist closer 
and closer to F and with increasing moment as the hole 
doping increases. For y = 0.3 the moment is nearly 0.5 
/i_B, except very close to F. A FM solution is thus being 
approached. 

Calculations for spirals with non-zero qz reveal much 
stronger dependence of the energy on the relative orien- 
tation of the Fe moments in adjacent FeAs layers than in 
LaFeAsO. This is due to the stronger dispersion of the 
As-pz hybridized bands which was discussed in Sec. IIIII 
In BaFe2As2, inter-layer nearest As neighbors are on top 
of each other and we find that the lowest-energy solu- 
tion is for nearest-neighbor layers having parallel AFM- 
ordered stripes. The energy for FM ordering between 
parallel stripes is 4 meV/Fe higher, and the energies for 
orthogonal stripes are intermediate. This is in accord 
with the experimental observations^ and results obtained 
from calculations for collinear spin arrangementsi^ 
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FIG. 6: (Color online) Dependences of the energy of 
LaFeAsOi_a;Fa; and Bai_2j/K2j,Fe2As2 on the angle a be- 
tween the moments of Fe nearest neighbors. The inset (turned 
45° with respect to the one in Fig. |4]| shows the spin arrange- 
ment in the Fe plane for a — 90° 



D. Applicability of the simple ji-j2 Heisenberg 
model 



The values of the stripe-ordered moment, M (X), the 
j2/ii ratio determined by the position of the local min- 
imum along MT, as well as the values of ji and j2 are 
collected in Table HI] for all electron and hole dopings con- 
sidered. Whereas M (X) and j2/ji exhibit clear trends 
with doping and position of the d^y band, ji and j2 scat- 
ter more. This enforces a conclusion that an incommen- 
surable SDW emerging upon electron doping is a typi- 
cal band-structure effect: The magnetic energies of the 
doped compounds can hardly be described by the simple 
ji-j2 Heisenberg model. 

In order to test the applicability of ji-j2 Heisenberg 
model for merely the parent compounds, we then per- 
formed calculations for spin structures with q = (1,0), 
but with the phases ipi for the two Fe sites in the tetrag- 
onal unit cell chosen in such a way that the angle, 
a ~ ipi — Lp2, between their magnetic moments varied 
from to 180°. As seen from the inset in Fig. [SI which 
is turned 45° with respect to the one in Fig. [H the Fe 
sites form two interpenetrating square sublattices. At 
q = (1,0), Fe moments in each sublattice are ordered 
antiferromagnetically as shown in the inset. The angles 
between the moment of an Fe and those on the two pairs 
of nearest neighbors, belonging to the other sublattice, 
are respectively a and tt — a and, hence, the Heisenberg 
energy is independent of a. The calculated dependencies 
of the energy on a are shown in Fig. [S) The spin ar- 
rangement corresponding to a=180° is exactly the same 
as the one generated by the spiral with q = (0, 1) and 
a=0. In other words, the spirals with a=0 and a=180° 
result in stripe order with the AFM chains parallel to x 
and y axis, respectively. The corresponding solutions for 



LaFeAsO are degenerate. In BaFe2As2 with the body- 
centered unit cell the variation in a is accompanied by a 
change of angle between the Fe moments in the adjacent 
FeAs layers. Due to the AFM alignment of the Fe mo- 
ments along the c axis, the solution with a=180° has a 
lower energy than that with a— having FM order along 
c. 

Our LSDA calculations show that for both compounds 
the energy depends strongly on the relative orientation of 
the Fe moments, with E{a=%Q) — E{a=Q) being compa- 
rable to the energy difference between the two collinear 
AFM solutions with q = (1,0) and q — (1, 1) discussed 
in connection with Fig. 21 The calculated results behave 
like E{a) — C siv? a, which does not appear in the ji-j2 
Heisenberg model but can be recovered by adding a bi- 
quadratic term proportional to (Si • Sj)^, with and Sj 
being spins on Fe nearest neighbors. 

The strong dependence of the energy on the relative 
orientation of the two AFM Fe sublattices points to the 
non-Heisenberg character of the interactions between Fe 
moments even in undoped LaFeAsO and BaFe2As2. The 
E{a) curves calculated for doped compounds (Fig. [S]) 
show a similar, although weaker, a dependence, espe- 
cially for Bao.6Ko.4Fe2As2. 

Thus, among the variety of spin configurations, which 
were degenerate in the ji-j2 Heisenberg model, a 
collinear spin arrangement with either q = (1,0) or 
q = (0, 1) is favored already at the level of the LSDA 
electronic structure. Such a magnetically ordered solu- 
tions lower the symmetry of the lattice from tetragonal to 
orthorhombic and lifts the degeneracy of Fe dyz and d^x 
states. This symmetry lowering is apparently responsible 
for the anisotropy of the exchange interactions calculated 
for stripe AFM order in Ref. US 



VI. NONINTERACTING SUSCEPTIBILITY 

With an Fe moment of '--'1.5 /ib, the exchange split- 
ting between the majority- and minority-spin Fe d states 
in the LSDA calculations is about 1.3 eV, which is as 
large as the dispersion of the paramagnetic Fe dxy and 
dyz and dxz bands over a significant part of the BZ, e.g., 
the dyz band in a tube around MY. Such a strong mag- 
netic perturbation dramatically changes the band struc- 
ture and the topology of the Fermi surface. The limit of 
a weak magnetic perturbation can be analyzed by study- 
ing the q and doping dependencies of the static, non- 
interacting linear-response susceptibility. Its imaginary 
part, Imxo(q) = Imxo(q, 0)/a;, is determined by 
the shape of and velocities on the Fermi surface (FS) 
and is a quantitative measure of FS nesting. The real 
part, Rexo(q) = Rexo(q, = 0), describes the response 
of the system to an infinitesimally small perturbation. 
In contrast to Imxo(q), the electronic states in a wide 
energy range around Ep may contribute to Rexo(ci)- 
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In the present work the noninteracting susceptibihty 

, , 1 ^ /„.(k + q)-^(k) 

XoK^.^) y 2^ e„,(k + q)-e„(k)+c^ + z<5 

X (k, nle-^i-nk + q, n') (k + q, n'|e'i"-|k, n) (1) 

was calculated using the linear response expressions given 
in Ref. l3J. Here, e„(k) is the energy of the nth band and 
/„(k) is the Fermi function. First, the imaginary part 
of xo{^.^^) was calculated in the 5^0 limit. Then, 
the real part was obtained by using Kramers-Kronig re- 
lations. The matrix elements (k + q, n'| exp(iqr)|k, n) of 
the perturbation were approximated by expanding the 
exponent inside each Fe sphere in Bessel functions and 
keeping only the spherically symmetric term proportional 
to jo (qr) . The expressions for the matrix elements were 
further simphfied by using jo(qr) ~ 1- Within this ap- 
proximation the contribution of each Fe sphere to the 
matrix element (k -I- q, n'| exp(iqr)|k, ti) is proportional 
to the overlap integral of the LMTO wave functions 5'k,n 
and ^k-(-q,n' inside the sphere. In other words, two states 
contribute to the susceptibility only if they have similar 
Fe partial-wave character. Although these approxima- 
tions are valid only at sufficiently small |q|, they do not 
affect the analysis of susceptibility peaks which may ap- 
pear due to the FS nesting. More details on calculation 
of yo(c|, u)) using the LMTO method can be found in 
Ref.lHi 

The contribution of a particular subset of electronic 
states to xo(q)W) can be discerned by calculating the 
matrix elements of the perturbation with all coefficients 
of the LMTO wave functions, except those which corre- 
spond to the chosen subset, set to zero. Due to the pres- 
ence of interference terms, such orbitally resolved contri- 
butions to the susceptibility are not additive. However, 
they allow us to discriminate those states which give the 
dominant contribution to xoi^.!^)- 

The susceptibilities of iron pnictides presented below 
were calculated starting from self-consistent non-spin- 
polarized electron densities and neglecting spin-orbit cou- 
pling. Thus, because of the degeneracy of Bloch states 
with different spin projections, xoi^.i^) calculated us- 
ing Eq. ^ describes the response of the system to a 
q-dependent spin as well as charge perturbation. A 
32x32x32 k mesh in the small tetragonal BZ was used 
in the calculations. 



A. LaFeAsOi_^F^ 

Nesting between the quasi-two-dimensional FS sheets 
in the iron pnictides were noted already in the earliest 
electronic-structure calculations Jii^i^S The FS cross sec- 
tions calculated for undoped LaFeAsO with and without 
the downwards shift of the Fe dxy energy were shown in 
Fig. [T] and their nesting pointed out in the caption to this 
figure, as well as in Sec. IIIIl 



The imaginary and real parts of the bare static suscep- 
tibility calculated for LaFeAsOi_a;Fa; with downshifted 
dxy energy are shown in the left-hand panel of Fig. [T] 
We first discuss the results for the undoped [x = 0) com- 
pound. In agreement with previous results"'^^'"'^^ the max- 
imum of Rexo(q) is found at q = X (1, 0). A peak of the 
imaginary part is, however, shifted away from X toward 
M. Analysis of the partial- wave resolved contributions 
to the susceptibility shows that the main contribution to 
Imxo(X) comes from nesting of M-centered dyz_a:z-like 
hole sheet (blue in the left-hand side of Fig. [T]) with 
the Y-centered electron sheet (green). This also gives 
the dominant contribution to the maximum of Rexo (X) . 
But also the dxy states (red) contribute significantly to 
Rexo(X) although their contribution to Imxo(X) nearly 
vanishes due to the ineffective nesting of the innermost 
F-centered hole sheet with the X-centered electron sheet. 
As q moves along XM, the hole and electron sheets start 
to touch when q k, (1, 0.13) and a peak of the dxy contri- 
bution to Imxo(q) appears at this nesting vector. This is 
responsible for the maximum of Imxo(q) along XM. The 
dxy contribution to Rexo (q) reaches its maximum at the 
same q. Since the weight of the Fe d^^'^-i and dx^-yi 
states in the bands crossing Ep is very small, they do 
not contribute to Imxo(q)5 whereas their contribution to 
Rexo(q) is almost constant in the whole q range. 

With electron doping, the hole sheets, centered at f 
and M, shrink and the electron sheets, centered at X 
and Y, grow. The hole and electron sheets no longer 
nest for q =X so that the susceptibility at X decreases 
rapidly with doping. Instead, peaks develop in both the 
imaginary and real parts of Xo(q) for q along XM for 
which the electron and hole sheets touch. This shift of 
the susceptibility peaks with increasing x correlates with 
the shift of the minimum of the E{q) curves calculated 
for spin spirals (Fig. 2]). The shift of the Rexo(q) peak 
along the XF line, accompanied by strong suppression at 
X, was noted in Ref. [ll|. 

Comparison of xo (q) calculated with (Fig. [7]) and with- 
out (Fig. [5]) shifting the Fe dxy energy shows that the 
susceptibility of undoped LaFeAsO is sensitive to the 
changes of the FS shown in Fig. [T] caused by the shift. 
Without the shift, the f-centered da;^-like hole sheet is 
larger and nests almost perfectly the X-centered electron 
sheet when q = X. This leads to Imxo(X) being the top 
of a sharp peak with the dominant contribution coming 
from the dxy states. The dxy contribution to Rexo(q) 
also has a sharp maximum at X. Since even the largest 
of the two M-centered dyz,zx hole sheets is smaller than 
the Y-centered electron sheet, nesting of these sheets does 
not contribute to Imxo(X). Their contribution increases 
as q shifts away from X, but remains weaker than the dxy 
one. Nevertheless, the dyz,zx contribution to Rexo(q) in 
the vicinity of the X point is comparable to the dxy one. 
The q dependence of the imaginary part of susceptibility 
is strongly affected by the change of FS nesting caused 
by the shift of the dxy bands, but the behavior of the real 
part is much more robust. Its maximum does not move 
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FIG. 7: (Color online) Imaginary and real parts of the bare static susceptibility for LaFeAsOi-i^Fa; and Bai_2yK2j/Fe2As2. 



away from the X point but looses its sharpness as the FS 
nesting becomes less perfect. 

The doping dependence of the susceptibility is only 
weakly affected by the shift of the dxy states. This may 
be one reason why different layered iron arsenides exhibit 
similar properties, in spite of the variation in their band 
structures. 



B. Bai_2i/K2;,Fe2As2 

As explained in Sec. IIIIl the FS of BaFe2As2 is simi- 
lar to that of LaFeAsO but all cylinders, except the F- 
centered d^y-like hole cylinder, are far more warped in 
kz direction. Besides, the bet symmetry mixes disper- 
sion into the {kx,ky) dispersions. Although the nesting 
for g = X is good when k^ = 7r/2c, it deteriorates at 
other kz- Nevertheless, a peak of Imxo(q), mostly due 
to the dxy states, is still present at X (right-hand side of 
Fig. 21), and the behavior of Rexo(q) for y — is qual- 
itatively similar to that for LaFeAsO. This comparison 
leads to a conclusion that independently of the fine de- 
tails of FS nesting, the formation of commensurate stripe 
order in the undoped compounds is preferable also in the 
limit of a weak magnetism. 

Hole doping makes the hole sheets grow and the elec- 
tron sheets shrink. The nesting with q =X is, nonethe- 



less, destroyed just as effectively as by electron doping. 
As a consequence, the real and imaginary parts of the 
susceptibility at X are strongly suppressed. In contrast 
to LaFeAs0i_2:Fa;, the peak in Rexo(q) shifts toward F, 
like the energy minimum calculated for spin spirals in 
Bai_2j/K2j,Fe2As2 (Fig. |4]) for the highest doping. The 
appearance of small-g, small-moment spin spirals corre- 
lates with the flat nonvanishing behavior of Rexo(q) near 

r. 



VII. CONCLUSIONS 

In conclusion, our LSDA total-energy calculations for 
spin spirals in LaFeAsOi-ajF^, and Bai_2aK2j,Fe2As2 us- 
ing the experimental crystal structures confirm that in 
the undoped compounds the minimum of the total en- 
ergy is reached at the wave vector q = X(1,0), which 
corresponds to stripe AFM order. The stability of this 
solution is, however, strongly affected by doping. With 
electron doping (x) exceeding 0.1 in LaFeAsOi-^^Fa;, the 
minimum becomes shallow and shifts toward M(l,l) to 
an incommensurate wave vector. This destabilization 
of the commensurate collinear stripe order by electron 
doping is a band-structure effect but not sensitive to 
the details of Fermi surface nesting. Hole doping (y) 
in Bai_2j/K23^Fe2As2 roughly continues the trend calcu- 
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FIG. 8: (Color online) Imaginary and real parts of the bare 
susceptibility calculated for LaFeAsOi-iFj^ without shifting 
of the d^cy energy. 



lated for decreasing electron doping in LaFeAsOi-ajFj^, 
that is, y —x. The energy gain due to stripe for- 
mation decreases with hole doping, but the minimum 
stays at X for y < 0.25. In both compounds, the de- 
viation of the q dependence of the energy from that of 
the classical Heisenberg model with nearest- and next- 
nearest-neighbor interactions becomes more pronounced 
with doping. 

We found that even in the parent compounds the total 
energy varies strongly when two interpenetrating AFM 
sublattices formed by the Fe ions are rotated with respect 
to each other. The dependence of the energy on the angle 
between the Fe moments in the two sublattices cannot 
be described by the simple j'i~j2 Heisenberg model, but 
may be reproduced by a biquadratic term proportional 
to (Si • Sj)^, which favors coUinear stripe AFM order. 

Although the LSDA for the experimental crystal struc- 
tures gives a stripe moments around 1.5 /is/Fe and con- 
comitant eV-large exchange splittings of degenerate ej,k 
and Sj'^u+q bands, linear-response calculations of the real 
and imaginary parts of the static noninteracting suscep- 
tibility based on the paramagnetic LDA band structure 
give similar results as the charge- and spin-self-consistent 
spin-spiral calculations. 
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